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Abstract 

We compare Newtonian three-flavor multigroup Boltzmann (MGBT) and (Bruenn's) 
multigroup flux-limitcd diffusion (MGFLD) neutrino transport in postbounce core collapse 
supernova environments. We focus our study on quantities central to the postbounce neutrino 
heating mechanism for reviving the stalled shock. Stationary-state three-flavor neutrino 
distributions are developed in thermally and hydrodynamically frozen time slices obtained from 
core collapse and bounce simulations that implement Lagrangian hydrodynamics and MGFLD 
neutrino transport. We obtain distributions for time slices at 106 ms and 233 ms after core bounce 
for the core of a 15 Mq progenitor, and at 156 ms after core bounce for a 25 Mq progenitor. For 
both transport methods, the electron neutrino and antineutrino luminosities, RMS energies, and 
mean inverse flux factors, all of which enter the neutrino heating rates, are computed as functions 
of radius and compared. The net neutrino heating rates are also computed as functions of radius 
and compared. 

Notably, we find significant differences in neutrino luminosities and mean inverse flux factors 
between the two transport methods for both precollapse models and for all three time slices. In 
each case, the luminosities for each transport method begin to diverge above the neutrinospheres, 
where the MGBT luminosities become larger than their MGFLD counterparts, finally settling to 
a constant difference maintained to the edge of the core. We find that the mean inverse flux 
factors, which describe the degree of forward peaking in the neutrino radiation field, also differ 
significantly between the two transport methods, with MGBT providing more isotropic radiation 
fields in the gain region. 

Most important, we find, for a region above the gain radius, net heating rates for MGBT that 
are as much as ~ 2 times the corresponding MGFLD rates, and net cooling rates below the gain 
radius that are typically ~0.8 times the MGFLD rates. These differences stem from differences in 
the neutrino luminosities and mean inverse flux factors, which can be as much as 11% and 24%, 
respectively. They are greatest at earlier postbounce times for a given progenitor mass and, for a 
given postbounce time, greater for greater progenitor mass. We discuss the ramifications these 
new results have for the supernova mechanism. 

Subject headings: (stars:) supernovae: general 
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1. Introduction 

Ascertaining the core collapse supernova mechanism is a long-standing problem in astrophysics. 
The current paradigm begins with the collapse of a massive star's iron core and the generation 
of an outwardly propagating shock wave that results from core rebound. Because of nuclear 
dissociation and neutrino losses, the shock stagnates. This sets the stage for a shock reheating 
mechanism whereby neutrino energy deposition via electron neutrino and antincutrino absorption 
on nucleons behind the shock reenergizes it (Bethe & Wilson 1985; Wilson 1985). 

The shock reheating phase is essential to the supernova's success, but it is precisely this 
phase that is difficult to simulate realistically. During shock reheating, core electron neutrinos 
and antineutrinos are radiated from their respective neutrinospheres, and a small fraction of 
this radiated energy is absorbed in the exterior shocked mantle. The shock reheating depends 
sensitively on the electron neutrino and antineutrino luminosities, spectra (best characterized by 
the RMS energies), and angular distributions in the region behind the shock (e.g., see Burrows 
& Goshy 1993, Janka & Miiller 1996, Mezzacappa et al. 1998b). These, in turn, depend 
on the neutrino transport in the semitransparent region encompassing the neutrinospheres, 
necessitating a neutrino transport treatment that is able to transit accurately and seamlessly 
between neutrino-thick and neutrino-thin regions. 

Various neutrino transport approximations have been implemented in simulating core collapse 
supernovae. The most sophisticated approximation, which naturally has been used in detailed 
one-dimensional simulations, is multigroup flux-limited diffusion (MGFLD; e.g.. Bowers & Wilson 
1982, Bruenn 1985, Myra et al. 1987). MGFLD closes the neutrino radiation hydrodynamics 
hierarchy of equations at the level of the first moment (the neutrino flux) by imposing a 
relationship between the flux and the gradient of the neutrino energy density (the zeroth moment). 
For example, 

_ cKdU 

^ " l/X + \dUjdr\/W,' 

where A is the neutrino mean free path, and Uu and Fi, are the neutrino energy density and 
flux (Bruenn 1985). [Other forms for the flux-limiter A can be found in Bowers & Wilson 
(1982), Levermore & Pomraning (1981), and Myra et al (1987).] Whereas the limits A — >■ and 
A — > oo produce the correct diffusion and free streaming fluxes, it is in the critical intermediate 
regime where the MGFLD approximation is of unknown accuracy. Unfortunately, the quantities 
central to the postshock neutrino heating mentioned earlier are determined in this regime, and 
given the sensitivity of the neutrino heating to these quantities it becomes necessary to consider 
more accurate transport schemes. Moreover, in detailed one-dimensional simulations that have 
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implemented elaborate MGFLD neutrino transport (e.g., see Bruenn 1993, Wilson & Mayle 1993, 

and Swesty & Lattimer 1994), explosions were not obtained unless the neutrino heating was 
boosted by additional phenomena, such as convection. This leaves us with at least two possibilities 
to consider: (1) Failures to produce explosions in the absence of additional phenomena, such as 
convection, have resulted from inexact neutrino transport. (2) Additional phenomena may be 
essential in obtaining explosions. 

Option (1) requires further comment. All investigators agree convection will occur during 
the shock reheating, explosion initiation phase in core collapse supernovae (Herant, Benz, & 
Colgate 1992; Miller et al. 1993; Herant et al. 1994; Burrows et al. 1995, Janka & Miiller 
1996, Mezzacappa et al. 1998b). Therefore, strictly speaking, all investigators agree the flow 
will not be spherically symmetric. However, although convection will certainly occur, it may 
play no significant role in initiating the explosion. It is with this distinction in mind that option 
(1) need be considered. For example, significant neutrino-driven convection was seen in recent 
multidimensional simulations employing one-dimensional MGFLD neutrino transport (Mezzacappa 
et al. 1998b); however, the angle-averaged shock radius, among other quantities, did not differ 
significantly from its one-dimensional counterpart, and no explosion was obtained. 

Ultimately, any successful model of core collapse supernovae will have to reproduce observables 
that clearly do not originate from spherically symmetric explosions, the most obvious of which is 
neutron star kicks. At this point, whether or not these kicks are generated during the initiation 
of the explosion or shortly thereafter is an open question. Note in this regard that simulations 
that have invoked multidimensional effects such as convection to explain such kicks have had 
difficulty generating, for example, adequate kick velocities. Mechanisms invoking convection, or 
aspherical neutrino emission resulting from convection, have not been able to produce kicks in 
excess of about 300 km/s, which therefore cannot account for the fastest pulsars — for example, 
PSR 2224+65, which has a velocity around 800 km/s — (Janka k Miiller 1994; Burrows & 
Hayes 1995, 1996). Moreover, definitive predictions of neutron star kick velocities from aspherical 
supernovae will require three-dimensional simulations. (The aforementioned simulations were 
carried out in two dimensions.) Recent simulations of neutrino-driven convection in two and three 
dimensions demonstrate that, as expected, longer-wavelength modes break up in three dimensions, 
rendering the angle-averaged flow qualitatively much closer to spherically symmetric (Knerr et al. 
1998; see also Miiller 1993 and Janka and Miiller 1996). In light of this, it is difficult to sec how 
the already low values for neutron star kick velocities obtained by invoking convection and/or 
convection-induced anisotropic neutrino emission during the explosion itself could be enhanced 
when these same simulations are carried out in three dimensions. 

In this paper, we compare three-flavor multigroup Boltzmann neutrino transport (MGBT) and 
(Bruenn's) MGFLD in spherically symmetric, hydrostatic, thermally frozen, postbounce proflles, 
with an eye toward quantities central to the postbounce neutrino heating mechanism for reviving 
the stalled shock. In particular, for both transport methods, we compute and compare the 
neutrino luminosities, RMS energies, mean inverse flux factors, and net heating rates as functions 
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of radius, time, and precollapse model. We then discuss the ramifications our results have for the 

supernova mechanism. This work is a continuation and extension of core collapse simulations 
(Mezzacappa and Matzner 1989, Mezzacappa and Brucnn 1993a, b,c), in which exact Boltzmann 
neutrino transport and multigroup flux-limited diffusion were compared. 

2. Initial Models, Codes, and Methodology 

We begin with 15 Mq and 25 Mq precollapse models S15s7b and S25s7b provided by Woosley 
(1995). The initial models were evolved through core collapse and bounce using one-dimensional 
Lagrangian hydrodynamics and MGFLD neutrino transport coupled to the Lattimer-Swesty 
equation of state (Lattimer & Swcsty 1991). The data at 106 ms and 233 ms after bounce 
for S15s7b and 156 ms after bounce for S25s7b were thermally and hydro dynamically frozen. 
Stationary-state neutrino distributions were computed for these profiles using both MGBT and 
MGFLD. 

The MGBT simulations were performed using BOLTZTRAN: a Newtonian gravity, 0{v/c), 
three-flavor, Boltzmann neutrino transport code developed for the supernova problem and used 
thus far for studies of stellar core collapse (Mezzacappa & Matzner 1989, Mezzacappa Sz Bruenn 
1993abc). The MGFLD simulations were performed using MGFLD-TRANS: a Newtonian gravity, 
0{v/c), three-flavor, MGFLD neutrino transport code, which has been used for both core collapse 
and postbounce evolution (e.g., Bruenn 1985, 1993). 

The MGBT simulations used 110 nonuniform spatial zones spanning radii from the origin to 
4744 km and 4673 km for model S15s7b at tpb = 106 ms and 233 ms respectively, and to 2096 
km for model S25s7b at tpt, = 156 ms. Twelve energy zones spanning a range between 5 and 225 
MeV were used to resolve the neutrino spectra. The MGFLD used the same spatial and energy 
grids. Simulations with 20 energy zones spanning the same energy range were performed with 
BOLTZTRAN; no changes in the results presented here were seen. 

For the MGBT simulations there is an added dimension: neutrino direction cosine. Because 
MGBT computes the neutrino distributions as a function of direction cosine and energy for each 
spatial zone, the isotropy of the neutrino radiation field as a function of radius and neutrino energy 
is computed from first principles. This is one of the key features distinguishing MGBT and MGFLD. 
Because the isotropy of the neutrino radiation field is critical to the shock reheating/revival, four 
Gaussian quadrature sets (2-, 4-, 6-, and 8-point) were implemented in the MGBT simulations to 
ensure numerical convergence of the results. 
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3. Results 

For electron neutrino and antineutrino absorption on neutrons and protons, the neutrino 
heating rate (in MeV/nucleon) in the region between the neutrinospheres and the shock can be 
written as 



where: Xn,p are the neutron and proton fractions; Aq = Ag = G\p{gy + 3g\)/TT{hc)^mB] 
Gp/i^c)^ = 1.166 X 10~^ GeV~^ is the Fermi coupling constant; p is the matter density; gy = 1-0, 
gA = 1.23; niB is the baryon mass; and L^^^p^, < E^^ p^ >, and F, F are the electron neutrino and 
antineutrino luminosities, RMS energies, and mean inverse flux factors, defined by 

27rc C 

Lv, = ^^'^'^Jf^ J dE^odpu^Kf^uJ, (4) 
fj,2._ !dE,Jp,^Elj 



JdE^Jp^^Elf 
,1 JdE,Jp,^Elf cU,, 



(6) 



y<2 



In equations (@)^®) / is the electron neutrino distribution function, which is a function of the 
electron neutrino direction cosine, Hu^, and energy, Ey^. In equation @, Uy^ and Fy^ are the 
electron neutrino energy density and flux. Corresponding quantities can be defined for the electron 
antineutrinos. Success or failure to generate explosions via neutrino reheating must ultimately rest 
on the three quantities defined in equations Both the MGBT and the MGFLD stationary 

state distributions were computed in the same thermally and hydro dynamically frozen matter 
configuration. 

In Figure 1, at 233 ms after bounce for model S15s7b, we plot the electron neutrino and 
antineutrino RMS energies, luminosities, and mean inverse flux factors as functions of radius for 
our (8-point Gaussian quadrature) MGBT and MGFLD runs. Energy-averaged electron neutrino- 
and antineutrino-spheres were located by calculating an energy-integrated neutrino depth defined 
by 

__ i:,dr' J,^dEydf^yE^yf/X 

J^dEyd^iyElf ' ^ ^ 
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and determining the radius at which r = 2/3. The neutrinospheres (at 57 km and 48 km, for 
electron neutrinos and antineutrinos, respectively), and the location of the shock (at 191 km), 
are indicated by arrows. The gain radius (neutrino-energy integrated), located at 98 km, is also 
marked by an arrow. For the electron neutrinos, the differences in RMS energies between MGBT 
and MGFLD are at most 2% throughout most of the region plotted, although MGBT consistently 
gives higher energies. The differences between MGBT and MGFLD antineutrino RMS energies are 
smaller, and neither transport scheme yields consistently higher values. It should be noted that 
we expect larger differences when a fully hydrodynamic simulation is carried out, with MGBT 
giving harder spectra (Mezzacappa and Bruenn 1993a, c; see also Burrows 1998). In a static 
matter configuration, differences that result from different treatments of the neutrino energy shift 
measured by comoving observers do not occur. 

Significant differences between MGBT and MGFLD are evident when comparing the neutrino 
and antineutrino luminosities and mean inverse flux factors. The luminosity curves for both 
electron neutrinos and antineutrinos coincide for both transport methods until the neutrinospheres 
are approached from below. Just below the neutrinospheres, the MGBT luminosities diverge upward 
from the MGFLD luminosities, differing by 7% (4% for antineutrinos) at the neutrinospheres. The 
root cause of this difference is that the MGFLD interpolation underestimates the flux in this region. 
After a decline from this maximum difference, the fractional difference grows from approximately 
3% at the base of the gain region to a constant difference of 6% beyond about 170 km. Similar 
behavior is exhibited by the antineutrino luminosities, with the same fractional differences, 3% 
and 6%, seen at the base of the gain region and near the shock, respectively. 

For the electron neutrinos, the fractional difference between < 1/F >mgfld and < 1/F >mgbt 
is 2%, 8%, and 12% at the neutrinosphere, gain radius, and shock, respectively. Just above the 
shock, the difference converges to 10%, and is maintained to the edge of the core. Focusing on 
the semitransparent region, < 1/F >mgfld is greater until the gain radius is approached from 
below; i.e., the MGFLD neutrino radiation field is more isotropic than the MGBT radiation field 
below these radii. At 80 km, as the gain radius is approached, MGFLD computes a sharp decrease 
in < 1/F >. Looking at Figure 2, where we plot the density and the sum of the MGBT electron 
neutrino and antineutrino luminosities for this time slice as functions of radius, it is evident as 
the gain radius is approached from below that the luminosity sum begins to turn over, marking 
the enclosure of the neutrino and antineutrino source. Therefore, for MGFLD, the accelerated 
transition to free streaming occurs not at the neutrinospheres, as might have been anticipated, but 
at a radius within which most of the neutrino and antineutrino source is enclosed. For example, 
for this time slice the electron neutrino luminosity at the neutrinosphere is only 76% of its peak 
value. 

In Figure 1, there is a second dip and a small spike in < 1/F >mgfld at 106 km and 163 km, 
respectively; < 1/F >mgbt is smooth through these radii. Again, examining Figure 2, the density 

shows discontinuities at 106 km and 163 km, which produce these features. The density profile 
flattens at 106 km and then drops precipitously. There is a corresponding flattening and sharp 
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drop in < 1/F >mgfld in this region. The density actually increases at 160 km, then immediately 
falls off. This results in an increase in < 1/F >mgfld at that radius, followed by a sharp decrease. 
In both cases, the isotropy of the MGFLD neutrino radiation field is altered by local conditions. 

For the electron antineutrinos, the same features are seen in < 1/F >mgfld- The fractional 
difference is 0%, 11%, and 11% at the neutrinosphcre, gain radius, and shock, respectively. The 
initial sharp decrease in < 1/F >mgfld occurs at a smaller radius. The antineutrino luminosity 
maximum, i.e., the point at which the antineutrino source is enclosed, is at a smaller radius. 

In Figure 3, wc plot the same three quantities for the earlier time slice in our 15 M0 model 
(tpb = 106 ms). The differences between MGBT and MGFLD are similar to those seen at the 
later postbounce time (tpb = 233 ms). The electron neutrino RMS energies are slightly higher 
for MGBT, again by about 2%. The differences in antineutrino RMS energies are again variable 
in both sign and magnitude, but are never more than 2%. The difference in luminosity is 11% 
(8% for antineutrinos) at the neutrinosphcre, 11% (7% for antineutrinos) at the gain radius, and 
settles to a constant difference of 9% (6% for antineutrinos) above 170 km. As in the later time 
slice, MGFLD underestimates the flux beginning below the neutrinospheres, extending everywhere 
above the neutrinospheres, which in turn results in a lower luminosity. 

For electron neutrinos, the fractional difference between < 1/F >mgfld and < 1/F >mgbt 
is 4%, 2%, and 17% at the neutrinosphere, gain radius, and shock, respectively. Most important, 
the same abrupt decrease in < 1/F >mgfld beginning just above the gain radius is evident. 
In Figure 4, we plot the density and the sum of the MGBT electron neutrino and antineutrino 
luminosities for tpb = 106 ms. The decrease in < 1/F >mgfld again occurs near the radius where 
the luminosity sum turns over: 128 km. There is also a small dip in < 1/F >mgfld at 172 km. 
Looking at Figure 4, the only significant dip in density occurs at 172 km, causing a local decrease 
in the isotropy of the MGFLD radiation field. 

For the electron antineutrinos, the difference between < 1/F >mgfld and < 1/F >mgbt 
is 3%, 2%, and 16% at the neutrinosphere, gain radius, and shock, respectively. The point at 
which < 1/F >MGFLD drops below < 1/F >mgbt is translated inwards, as expected, towards the 
antineutrino luminosity maximum. For both electron neutrinos and antineutrinos, a difference of 
about 13% is maintained with increasing radius above the shock. 

Figure 5 contains the same information as Figures 1 and 3, but for our 25 Mq model at 156 ms 
after bounce. The similarities between the 25 Mq and the 15 Mq results are striking, considering 
the marked difference in core structure. The identical trend in RMS energies is again seen: there 
are small differences for the electron neutrinos, but consistently higher values are obtained with 
MGBT (again, ~ 2% higher); there are smaller, sometimes oscillating, differences between MGBT 
and MGFLD for the electron antineutrinos. The absolute value of the neutrino luminosities at 
the gain radius is almost a factor of 2 greater than in either of the 15 Mq slices, yet the relative 
differences between the luminosities computed by the two transport methods are similar. The 
electron neutrino luminosities differ by 9% (5% for antineutrinos) at the neutrinosphere, 11% (9% 
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for antineutrinos) at the gain radius, and 9% (7% for antineutrinos) at the shock. The difference 
at the shock is maintained to the edge of the core. Again, this difference in luminosity is caused 
by an underestimation of the flux by MGFLD. 

Among the three slices considered, the 25 Mq core gives rise to the most dramatic differences 
in < 1/F >. For electron neutrinos, we find fractional differences of 1%, 2%, and 24% in 

< 1/F > at the neutrinosphere, gain radius, and shock, respectively. There is a sharp decrease in 

< 1/F >MGFLD at 109 km. As in the 15 Mq case, this change in < 1/F >mgfld is correlated with 
the enclosure of the neutrino source. This correlation is evident in Figure 6, where we plot the 
density and sum of the MGBT electron neutrino and antineutrino luminosities for this time slice. 
Note that the MGBT luminosity sum begins to turn over near 109 km, coincident with the sharp 
decrease in < 1/F >mgfld- 

There are other precipitous drops in < 1/F >mgfld at 125 km and 162 km. Also apparent in 
Figure 6 are precipitous drops in density at these radii. These drops in density produce changes 
in the local neutrino radiation field computed by MGFLD. 

For electron antineutrinos, < 1/F >mgfld exhibits similar structure, with fractional 
differences of 0%, 1%, and 19% at the neutrinosphere, gain radius, and shock, respectively. 

< 1/F >MGFLD for antineutrinos also contains three dips: < 1/F >mgfld drops below 

< 1/F >MGBT at 105 km, where most of the antineutrino source is enclosed, and additional dips at 
125 km and 162 km are also evident, again resulting from the density dips visible in Figure 6. For 
both electron neutrinos and antineutrinos, constant differences ~ 12% are seen above the shock. 

Because each of the quantities plotted in Figures 1, 3, and 5 is consistently greater for 
MGBT (while this is not strictly true for the antineutrino RMS energies in our stationary state 
comparisons, in a fully dynamical simulation these energies will be consistently higher for MGBT 
[Mezzacappa and Bruenn 1993a,c; see also Burrows 1998]), and because the neutrino heating rate 
is proportional to each of them, MGBT yields a significantly higher heating rate. As an example, 
just above the gain radius for model S15s7b at tpb = 233 ms, at the point of peak heating, MGBT 
yields a heating rate from neutrino absorption that is (102%)^ x 110% x 112% of the MGFLD rate. 
Note in equation (^) that the heating rate depends linearly on both the neutrino luminosity and 
the mean inverse flux factor. A reliable determination of the heating rate in and around the gain 
region therefore depends on a realistic solution of the transport equation in which both of these 
quantities are determined self consistently. 

In Figure 7, for MGBT and MGFLD, we plot the net neutrino heating rates as functions 
of radius for model S15s7b at tpb = 233 ms. (As described in Section 2, the results from four 
Gaussian quadrature sets are plotted to demonstrate numerical convergence.) These rates include 
the contributions from both the electron neutrinos and antineutrinos, and were computed using 
the following formulae: 
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{de/dt)i = c I EldE,[^y\f^ - - ^"Dypihcf (8) 

where e is the internal energy per gram; ji are the electron neutrino or antineutrino 

energy, zeroth angular moment, absorption mean free path, and emissivity, respectively; i = 1 
corresponds to electron neutrinos, and i = 2 corresponds to electron antineutrinos. Only the 
contributions from neutrino emission and absorption were included: in our models, contributions 
from neutrino-electron scattering and other processes contribute only a few percent to the neutrino 
heating rate at and before 233 ms after bounce; they become more important (~15%) at later 
times. The (8-point Gaussian quadrature) MGBT simulation yields a net heating rate just above 
the gain radius that is ~1.3 times the MGFLD rate, and a net cooling rate below the gain radius 
that is consistently ~0.8 times the MGFLD rate. Comparing the net heating rate to the luminosity 
sum in Figure 2, the extent of the gain region (from 98 km to 180 km ) is well correlated with 
the region between the luminosity maximum and the radius where the luminosity levels off to a 
constant value (from 101 km to 170 km). Note also that the gain radius is located at a smaller 
radius for MGBT, and consequently, the size of the net heating region below the shock is somewhat 
larger. 

For (8-point Gaussian quadrature) MGBT and MGFLD, Figure 8 shows the net heating 
curves for S15s7b at an earlier postbounce time, tpb = 106 ms. The region between the 
neutrinospheres and the shock is a bit smaller: The shock is approximately at the same radius, 
but the energy-averaged electron neutrinosphere is closer to the shock (at 84 km, versus 58 km 
for tpb = 233 ms). Similar net heating profiles arc seen, but the differences between MGBT and 
MGFLD are even greater. MGBT yields a net heating rate that is ^2 times the MGFLD rate just 
above the base of the gain region. Again, the cooling rate is consistently '^0.8 times the MGFLD 
rate below the gain radius. The correlation between the gain region and the region between the 
luminosity-sum maximum and the point at which the luminosity sum is constant is also quite 
strong for this time slice. The gain radius (103 km) is just below the luminosity maximum (107 
km); also, the luminosity levels off at a radius (173 km) near the upper edge of the gain region 
(182 km). 

Figure 9 contains the net heating curves for our 25 Mq model at = 156 ms. Again, in this 
case the (8-point Gaussian quadrature) MGBT net heating rate is ~2 times the MGFLD rate just 
above the base of the gain region, and the net cooling rate below the gain radius is consistently 
~0.8 times the MGFLD rate. The gain radius (111 km) and the luminosity-sum maximum (108 
km) are strongly correlated, as are the radius marking the upper extent of the gain region (180 
km) and the radius at which the luminosity sum levels off (170 km). Although the differences in 
the net heating rate are comparable to those seen in the tpb = 106 ms slice for model S15s7b, they 
occur at a significantly later postbounce time. At a given postbounce time, MGBT provides the 
greatest enhancement in net heating for larger progenitor masses. 
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4. Summary, Discussion, and Conclusions 

Comparing three-flavor MGBT and three-flavor MGFLD in postbounce supernova 
environments, we find that MGBT leads to a significant increase/decrease in the net heating/cooling 
rate, particularly above/below the gain radius. The MGBT net heating rate can be as much as 
~2 times the MGFLD net heating rate above the gain radiiis, with net cooling rates that are 
typically ~0.8 times the MGFLD rate below the gain radius. These differences stem primarily 
from differences in the neutrino luminosities and mean inverse flux factors; the heating rate is 
linearly proportional to both these quantities, and differences in both add to produce a significant 
difference in the net heating rate. 

In Figure 10, we plot the sum of the electron neutrino and antineutrino luminosities computed 
in the MGFLD S15s7b dynamic simulation for several different postbounce times. It is important 
to note that the total luminosity changes by ^5-15% between 100 km and 200 km on time scales 
^30 ms. Moreover, the light crossing time between 100 km and 200 km is ~l/3 ms. Therefore, 

the neutrino source in our simulations changes on time scales that are two orders of magnitude 
greater than the time scales on which stationary state is established in this region. This suggests 
our stationary state results closely reflect what will occur in dynamic simulations. 

We also observe that the differences in the net heating rate are greatest at earlier postbounce 
times for a given progenitor mass, and greater at any given postbounce time for greater progenitor 
mass. This is illustrated in Table 1. The increase in net heating with increased progenitor mass is 
advantageous because of the slower fall-off in the preshock accretion ram pressure. 

These results have at least two important ramifications for the supernova mechanism: 

(1) With the dramatic increase in net heating above the gain radius, which is seen in all 
of our postbounce slices, it may be possible to obtain explosions in one dimension without 
multidimensional effects such as convection; this will be aided by the decrease in net cooling 
below the gain radius. It should be noted that our postbounce slices come from full radiation 
hydrodynamics simulations implementing MGFLD that marginally failed to produce explosions 
(Bruenn 1993). The marginality of Bruenn's simulations is an important motivating factor in 
comparing our MGBT results solely with his MGFLD results. All else being equal, increases in 
net heating of the magnitude documented here would most likely lead to explosions. However, 
simulations coupling MGBT and the core hydrodynamics must be carried out in order to compute 
any feedbacks. It remains to be seen whether the MGBT heating rate will remain sufficiently high 
to generate an explosion. Also, the effects of using general relativistic gravity, hydrodynamics, 
and neutrino transport must be explored, especially if explosions are obtained in the Newtonian 
limit. For example, the redshifted neutrino energies, the smaller gain region, and the greater 
infall velocities in the gain region will most likely conspire to make explosions more difficult to 
obtain (Mezzacappa et al. 1998b, DeNisco, Bruenn, & Mezzacappa 1998, Bruenn, DeNisco, & 
Mezzacappa 1998). 
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(2) With the dramatic increase in net heating occurring near the base of the gain region, 
we anticipate that MGBT coTipled to two-dimensional hydrodynamics will yield more vigorous 
neutrino-driven convection than was seen in Mezzacappa et al. (1998b), where two-dimensional 
hydrodynamics was coupled to one-dimensional MGFLD. The combination of increased net heating 
and more vigorous neutrino-driven convection would be more favorable for shock revival. 

In closing, our results are promising, and their ramifications for core collapse supernovae 
and, in particular, for the postbounce neutrino-heating, shock-revival mechanism await one- 
and two-dimensional dynamical simulations with MGBT coupled to the core hydrodynamics. 
One-dimensional simulations are currently underway, and we plan to report on them soon. 
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7. Figure Captions 

Fig. 1. — At = 233 ms, for model S15s7b, wc plot the electron neutrino and antineutrino RMS 
energies, luminosities, and mean inverse flux factors versus radius for both Boltzmann neutrino 
transport and Bruenn's MGFLD. The location of the energy-averaged electron neutrinosphere and 
anti-neutrinosphere and the locations of the gain radius and shock are indicated by arrows. 

Fig. 2. — The density and MGBT electron neutrino and antineutrino luminosity sum are plotted 
versus radius for model S15s7b at tpb = 233 ms. 

Fig. 3. — At tpb = 106 ms, for model S15s7b, we plot the electron neutrino and antineutrino RMS 
energies, luminosities, and mean inverse flux factors versus radius for both Boltzmann neutrino 
transport and Bruenn's MGFLD. The location of the energy-averaged electron neutrinosphere and 
anti-neutrinosphere and the location of the gain radius and shock are indicated by arrows. 

Fig. 4. — The density and MGBT electron neutrino and antineutrino luminosity sum are plotted 
versus radius for model S15s7b at tpb = 106 ms. 

Fig. 5. — At tpb = 156 ms, for model S25s7b, we plot the electron neutrino and antineutrino RMS 
energies, luminosities, and mean inverse flux factors versus radius for both Boltzmann neutrino 
transport and Bruenn's MGFLD. The location of the energy-averaged electron neutrinosphere and 
anti-neutrinosphere and the location of the gain radius and the shock are indicated by arrows. 

Fig. 6. — The density and MGBT electron neutrino and antineutrino luminosity sum are plotted 
versus radius for model S25s7b at tpb = 156 ms. 

Fig. 7. — At tpb = 233 ms, for model S15s7b, the net neutrino heating rates are plotted versus 
radius for both Boltzmann neutrino transport and MGFLD. The location of the energy-averaged 
electron neutrinosphere and the location of the shock are indicated by arrows. The results from 
four Gaussian quadrature sets are plotted to demonstrate numerical convergence. 

Fig. 8. — At ipb = 106 ms, for model S15s7b, the net neutrino heating rates for MGFLD and (8-point 
Gaussian quadrature) Boltzmann transport are plotted versus radius, along with the locations of 
the energy-averaged electron neutrinosphere and the shock. 

Fig. 9. — At ipb = 156 ms, for model S25s7b, the net neutrino heating rates for MGFLD and (8-point 
Gaussian quadrature) Boltzmann transport are plotted versus radius, along with the locations of 
the energy-averaged electron neutrinosphere and the shock. 

Fig. 10. — The sum of the electron neutrino and antineutrino luminosities from the MGFLD dynamic 
run for model S15s7b are plotted versus radius. The luminosities from five different postbounce 
times are shown. 
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8. Tables 



Table 1: Maximum Net Heating/Cooling Rates 
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